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- _Init ialize random number seed 
nitialize random number seed 
BASIC Random Number Generator 
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TITLE Basan gyoon - BASIC Random Number Sysgert 
~IDENT /1-005/ ; File: BASRANDOM.MA 
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COPYRIGHT (c) 1978, 1980, 1982, 1984 B 
DIGITAL EQUIPMENT 1 CORPORATION, ° MAYNARD. MASSACHUSETTS. 
ALL RIGHTS RESERVED 


® 
® 
® 
® 
® 
THIS SOFTWARE I$ FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED  * 
ONLY IN JACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER * 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY * 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
TRANSFERRED. . 
® 
THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT ®* 
CORPORATION. * 
® 
® 
® 
® 
® 
* 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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FACILITY: BASIC Support Library 
ABSTRACT: 


This module provides iB he oe for all BASIC random number 
functions. The see in OWN storage in this module. 
There are 2 entry ALK, Phe first initializes the seed 
based on the time of day to prevent getting the same sequence 
each time a program is run. The second picks up the seed 

and generates a random number. 

This is a general conden number generator. It is 

of the multi sph teat’ ve congruential type, and hence 

is fast, although prone to certain classes of 

non-random sequences. 
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ENVIRONMENT: any access mode, normally user mode 
modular, non-ASf reentrant procedure 
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AUTHOR: R. Will, CREATION DATE: 29-Nov-78 
MODIFIED BY: 
R. Will, 29-Nov-78: VERSION 01 
nal crept ion, 
PSECT dire fives. JBS 22-DEC-78 
4-JUN-1979 


Add BASSERAND® INIT, JBS 
Don't cause ifteger BI wes 05-AUG-1979 
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1-005 = Start the generator at 1(€5983F7, so it 
first time it is called. JBS 61-SEP-1 


gogsn’ t give 0 the 
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BASSRANDOM - BASIC Random Number Support ot 14 -1984 AX/VMS 04-00 P 
1-005 DECLARATIONS ad Breer bbs 90595599 PBASK TS. Cae TBALGANSOM.mans? 2% (3) 
6¢ -SBTTL DECLARATIONS 
§ 64 ; 
65 ; INCLUDE FILES: 
00 66 : 
8 of : NONE 
09 $9 ; EXTERNAL DECLARATIONS: 
000 71 -DSABL GBL 3; Prevents undeclared symbols 
000 Me i from being automatically global 
4 i? : -EXTRN SYSSGETTIM 3; System services, current time 
000 75 ; MACROS: 
0000 6 3 
0000 73 NONE 
0000 78 ; 
0000 79 ; EQUATED SYMBOLS: 
0000 80 ; 
0000 81; 
0000 Ht 3; OWN STORAGE: 
0000 H ; 
0000 tH Be -PSECT _BASSDATA,PIC,NOSHR,WRT,LONG,NOEXE 
1C05983F7 Babe Hy SEED: -LONG ‘*X1C5983F7 ; Initial value of seed 
0006 88: 
Bae + 3; PSECT DECLARATIONS: 
00000000 91 : -PSECT _BASSCODE,PIC,USR,CON,REL,LCL,RD,SHR,NOWRT,LONG,EXE 
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BASSRANDOM = BASIC Random Number Support 16-SEP-1984 00:01: AX/VMS Macro V04-00 Page 
1 5 BASSRANDOMIZE = Initialize random numbe eee 98:91:90 EBASRTL. SRE IBASRANDOM. MAR: 1 ° Ry 
0 $4 “ ~SBTTL BASSRANDOMIZE - Initialize random number seed 
09 32 3 FUNCTIONAL DESCRIPTION: 
§ 39 ; This routine initializes ihe seed by getting the current system 
Bp 4 ; time and adding the low 32 bits to the high 32 bits to get the seed. 
000 100 ; CALLING SEQUENCE: 
0009 18) : CALL BASSRANDOMIZE ( ) 
0009 198 : INPUT PARAMETERS 
8009 188 3 
0000 106; NONE 
0000 107; 
0000 108 ; IMPLICIT INPUTS: 
0000 109; 
0000 110; NONE 
0000 111; 
0000 We 3; OUTPUT PARAMETERS: 
0000 113; 
0000 114; NONE 
0000 115; 
0000 116 ; IMPLICIT OUTPUTS: 
0000 117 ; 
0000 i18; SEED.wlu 
0000 119; 
0000 120 ; FUNCTION VALUE: 
0000 15) 3; COMPLETION CODES: 
0000 1 ; : 
0000 123; NONE 
B38 124 ; 
000 125 ; SIDE EFFECTS: 
0000 126; : 
0000 127; SEED is altered 
0000 128 ; 
0000 \¢? 3e° 
0000 130 
0000 bo03 13 -ENTRY BASSRANDOMIZE , “M<> ; Entry point 
7E 7C 0002 138 CLRQ -(SP) ; space on stack for time 
0004 134 S$GETTIM_S (SP) ; get current time 
000D 135 3 as top 2 words on stack 
0O0000000°EF 8—& 8E& C1 00D 136 ADDL3 (SP)+,(SP)+,SEED ; add the two words o: 
0015 137 3; current time to get 
0015 138 3; randomized seed 
0015 139 ; and clean up stack 
04 0015 140 RET 
| 
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-SBTTL BASSSRAND_INIT = Initialize random number seed 
FUNCTIONAL DESCRIPTION: | 
This routine sets the seed to 1(5983F7. It is used by the RUN 
l 


2 
om 


command when the environment is initialized. 
CALLING SEQUENCE: 
CALL BASSSRAND_INIT ( ) 
INPUT PARAMETERS: 
NONE 
IMPLICIT INPUTS: | 
NONE 
OUTPUT PARAMETERS: 
NONE 
IMPLICIT OUTPUTS: 
SEED.wlu 


FUNCTION VALUE: 
COMPLETION CODES: 


NONE 
SIDE EFFECTS: 
SEED is altered 


Be Se Se Se Se Se Ge Se Se Se Se Se Ge Se Se Se Se Se Se Se Se Se Se Fe Sse Se Se Se Se Fe Ge Ss Fe Ge Se Se 
+ 


0000 
QOO00000"EF  1C5983F7 BF pe 


-ENTRY BASSSRAND_INIT , “M<> ; Entry point 


MOVL #°x1C5983F7, SEED : Initialize the seed 
RET 3 and return 
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Number Support 
BASIC Random Number Gener 


ae ~SBTTL BASSRND_F_R1 = BASIC Random Number Generator 
; FUNCTIONAL DESCRIPTION: 
; RANDOM = BASIC Pseudo Random Number Generator 
Call the function BASSRND_F ( ) to obtain the 
next pseudo-ranaom number. The seed is updated by 
the function automatically as a side effect. The 
distributed in the range 


There are no coete tes tons on the seed, although 


it should be initialized to different values on 


separate runs. 


work: 


16-SEP-19 
6-SEP=19 


to multiply by the number 
if a uniform choice among 


(New 


Be 90:95:71 


AX/VMS Macro v04-00 
BASRTL.SRCJBASRANDOM.MAR; 1 


result is a floating point number that is uniformly 
0.0 inclusive to 1.0 exclusive. 


The algorithm used is to update the seed as: 
seed = 69069. * seed + 1 
and then to convert the seed to floating point. 


(mod 2**32) 


Note, because the result is never 1.0, a simple 
to get a uniform random integer selector is 


f cases. For example 
situations is to be 


Encyclopedia of Computer Science 
edited by Anthon 
Petrocelli 
; CALLING SEQUENCE: 
JSB result.wf.v = BASSRND_F_R1 ( ) 
; INPUT PARAMCTERS: 


NONE 


Ralston 
ork, 1976) 


made, then the following BASIC statement will 
ON 14INT(5.O*RND) 1,2,3.4,5 


Note that the explicit INT is necessary before 
in order to avoid a possible meaynetae 
during the normalization after the floating add. 


This is a 
of the multi 
is fast, alt 
non-random sequences. 


eneral random number generator. It is 
icative congruential type, and hence 
ugh prone to certain classes of 

This non-random behavior 
typically arises when considering triples of 
numbers generated by this method. 


For more information on congruential generators, 


Random Number Generation (pp. 1192-1197) 
by G. Marsaglia 
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| BASSRANDOM = BASIC Random Number Support 16-SEP-1984 00:01:20 \VAX/VMS Macro Vv04-00 
bet BASSRND_F_R1 = BASIC Random Number Gener ets 98:92:48 EBASRTL. SREIBA RANDOM .MAR; 1 
| ¢ 348; IMPLICIT INPUTS: 
| 4 45; 
4 44 ; SEED. mlu 
$ 4 45 ; 
4 246 ; OUTPUT PARAMETERS: 
| 024 47 ; 
024 48 ; NONE 
024 49 ; 
0024 50 ; IMPLICIT OUTPUTS: 
0024 2} ; 
0024 ¢ ° SEED.mlu 
8 4 29 3 
bp 2¢ 3 FUNCTION VALUE 
0024 56 ; Recurns in RO a single-precision floating point value between 
0024 of 3 0.0 inclusive and 1.0 exclusive. 
0024 58 ; 
002’ 59 ; SIDE EFFECTS: 
0024 60 ; 
0024 61 ; The value of SEED is altered 
Base $¢ 3 
0024 63 ;-- 
0024 64 BASSRND_F_R1:: ;JSB entry 
8 4 06 ; Do the multiply in a way that wilt not cause an integer 
4 r ef 3 overflow. 
01 Q0010DCD 8F epaeeeedh ** 7A + . 269 EMUL SEED ,#69069,41,R0 ; compute product in RO/R1 
00000000'EF 50 00 0038 4 MOVL RO,SEED ; ignore high-order bits of product 
003 72 ;+ 
One g ; The next instructions convert the seed from une igned integer 
0 3 74 ; to floating point in the range 0.0 to 1.0 exclusive. 
boss 57e* 
50 O0000000'EF 18 08 EF 0038 77 EXTZV #8,#24,SEED,RO :Get the most significant bits 
0041 re ; of the seed in the range 
0041 7 3.0 1. (2"#24)-1 
50 50 4E€ 0041 0 CVTLF =RO,RO :Convert to floating without 
0044 1 3; rounding. The result is 
0044 ; 3 positive and in the range 
0044 3 0.0 .. (2.0"*24)-1.0 
ape 4 
044 5 3+ 
0044 § 3; If this were to be placed as an inline expansion, then 
apoE i 3: MULF #*x00003480,R0 could replace the next two instructions. 
bee $89 
05 13 0044 90 BEQL 10$ :If zero, already correct 
50 0C00 BF A2 0046 291 SUBW = #24€7,RO :DIVF #*F2.0%824 
048 35 ; the result is now in the 
Ree 37 ; range 0.0 .. 1.0 exclusive 
05 0048 95 10$: RSB 
04C 38 
04C 9 ~ END 
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BASSRANDOM 
Symbol table 


BASSSRAND INIT 0000 16 RG 

BASSRANOORIZE it 0 RG 

BASSRND_F_R 000024 RG =O 
000 OR 01 
eererene 00 


SEED 
SYSSGETTIM 
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! Psect synopsis ! 
ree oes ee maa eam + 


PSECT name Allocation PSECT No. Attributes 

- ABS. 00000000 < 0.) 00 ¢ QO.) NOPIC US 

_BASSDATA 00000004 < 4.) O1¢ 1.) PIC US 

_BASSCODE Q0000004C <¢ =e ce PIC US 
fewer amr ee swe e sn wee wmaace + 
H Performance inuicators ! 

Phase Page faults CPU Time Elapsed Time 

Initialization — 3 80:00: 00.08 00:00:00.65 

Command processing 11 00:00:00.4 00:00:03.67 

Pass 1 102 BF 8: Bn “Be Bo 80 06 “RB 

Symbol table sort 0 00:00:00.00 00:00:00.00 

Pass 2 61 00:00:00.56 00:00:01.19 

Symbol table output ¢ 00:00:00.01 00:00:00.01 

Psect synopsis output 00:00:00.02 00:00:00.02 

Cross-reference output 9 00:00:00.00 00:00:00.00 

Assembler run totals 31 00:00:01.76 00:00:08.22 


The working set Limit was 1050 pages. 


3063 bytes (6 pages) of virtual memory were used to buffer the intermediate code. 
non-local ane 1 local symbols. 


There were 10 pages of symbol table space allocated to hold 
97 source Lines were read in Pass 1, produc ing 16 object records in Pass 
m 


peges of virtual memory were used to define acros. 
FTO me ON OR Om mM OR ee + 
! Macro Library statistics : 
mer or tO Bee BOE OO MO ee 


Macro library name Macros defined 
-$255$DUA28: CSYSLIBISTARLET.MLB; 2 2 
12 GETS were required to define 2 macros. 


There were no errors, warnings or information messages. 


RECIMEE W8E 


R CON 
R CON 
R CON 


SS 
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SHR = EXE 
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SHR NOEXE NORD 
HR NOEXE RD 


v04-00 p 
a 


NOWRT NOVEC BYTE 
WRT NOVEC LONG 
RD NOWRT NOVEC LONG 


MACRO/ENABLE=SUPPRESS10N/DI SABLE=(GLOBAL , TRACEBACK) /LIS=L1S$:BASRANDOM/OBJ=0BJ$:BASRANDOM MSRC$:BASRANDOM/UPDATE = (ENHS: BASRANDOM) 
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